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Abstract 

Damping on an object generally depends on its conformation (shape size etc.). We consider 
the Langevin dynamics of a model system with a conformation dependent damping and generalize 
the fluctuation dissipation relation to fit in such a situation. We derive equilibrium distribution 
function for such a case which converges to the standard Boltzmann form at the limit of uniform 
damping. The results can have implications, in general, for barrier overcoming processes where 
standard Boltzmann statistics is slow. 
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INTRODUCTION 



The Langevin dynamics for equilibrium processes in general includes a uniform damp- 
ing constant T over space for a particle that bears a relationship with the strength of the 
Gaussian white noise given by the expression yfZFKBT ', where, Kb is Boltzmann-constant 
and T is the temperature of the bath with which the particle is in equilibrium. In such a 
modelling, the average stochastic force experienced by the particle under consideration has 
to be a function (as mentioned above) of the damping constant of the particle to maintain 
an energy balance over mesoscopic time scales in equilibrium. The other important thing 
is that the T is generally structure dependent. For example, two spheres of different stokes 
radii in equilibrium with the same heat bath will have two different T and would see the same 
temperature T of the bath. Keeping this in mind, an object which can have conformation 
fluctuations in contact with a heat bath can in general have a damping constant dependent 
on its conformational states. In the following, we would do an analysis for this kind of a 
system with a generalization of the fluctuation-dissipation (FD) relation in order to ensure 
a stationary probability distribution in the conformational space of the system. This is a 
necessity to enforce no current condition in the particularly inhomogeneous conformational 
space of the system. The equilibrium distribution that we get is a generalization of Boltz- 
mann distribution function which becomes standard Boltzmann distribution at the limit of 
constant damping or on an average over the whole conformation-space. 



MODEL 

We consider a dimer model somewhat similar to the one we considered in Ref. to show 
a special case of having directed transport under equilibrium conditions. In fact, the directed 
transport under equilibrium conditions would also have a general proof as a byproduct of 
our present exercise. The one-dimensional model in the over-damped regime reads as 



r 2 (*i - X2 )x 2 = *A + r 2 ( Xl - ^ <T ^ X2)> vM (2) 

In the above expressions, Z = x\ —x 2 is the coordinate that captures various conformations of 
the system, U(Z) is the internal field that keeps particles talking (e.g. keeps them bound and 



takes care of the excluded volume interaction) and there is no external global field present. 
The noises have been taken explicitly in accordance with the generalized FD relation we 
talked about where < Y^Z) > is the spatial average of the damping constant which can be 
evaluated in a self consistent manner. Note that, for a uniform V we will recover the standard 
FD relation and the form of the standard FD relation is also there globally in terms of the 
average Y(Z). rjiS are the Gaussian white noise of unit strength. Our subsequent analysis 
would justify the reason for such a generalization of the FD relation. 

Moving to the conformation and centre of mass (X) coordinates we can rewrite the above 
model as 

._ Y 1 {z) + v 2 {z)dU{z) f^r^ ( Vl (t) m (t) \ 

z " " W) dz + V^bT ^ <ri(z)> " ^ <V2{Z)> ) ® 

■ Y l {Z)-Y 2 {Z)dU{Z) [KJT ( 7/! (t) y 2 (t) \ 

2Y 1 (Z)Y 2 (Z) dZ + V 2 r l( Z)> %/< Y 2 {Z)>) U 

Note that, the second moment of the noise terms (£z and £x) in the Eq.3 and 4 respectively 
are 

. c (t V (t W 2K B T(< Ti(Z)> + <Y 2 (Z) >) 

< >= <Vi{z) ><T2{Z)> K-t 2 (5) 

< P (t V (t ) >- KbT{< Tl{Z) > + < T2{Z) >] 5 (6) 

< M*0Mfc) >- 2 < ri(z) >< r 2( z) > — 5t ^ (6) 

where, the first moments vanish following the properties of rjiS. There is also a cross correla- 
tion between £ z and £x which we do not require to consider here because of the decoupling 
of the Z dynamics from that of X. Let us consider the short hand notations 

_ Y 1 (Z) + Y 2 (Z) 

u r 1 (z)r 2 (z) 

and 

_ Y 1 (Z)-Y 2 (Z) 
^ X 2Y l (Z)Y 2 (Z) 

and F(Z) = 

Now, the system to have a stationary probability distribution P(Z) that makes the aver- 
age current in the (inhomogeneous) conformation space < Z > vanish, we get the condition 
< z >=< (zF(Z) >= 0. This means, 

< Z >= f dZ( z F(Z)P(Z) = 0, (7) 

J a 



where the integration limits are at the two extreme zeros of the P(Z) beyond which the 
conformations are not visited. Let us impose a sufficient condition for this integral to vanish 
namely 

dP{Z) 



d{Z) 

which would actually solve to give us 



C( Z F(Z)P(Z), (8) 



P{Z) = iVexp (C J dZCzF(Z)), (9) 

where C is a constant that would make the exponent dimensionless and we would find it 
out from the corresponding Fokker-Planck (FP) equation. 



FOKKER-PLANCK EQUATION 

The FP equation for the Langevin dynamics of the conformation coordinate i.e. 

Z = CzF(Z) + £ z (t) (10) 
can be found out following the standard procedure mentioned in j^j], using the equation 

^ = -*<*»*>• ("> 

where the <> indicates a noise average for the noise distribution -P[£z(t)] = exp (— J t */ dt J^j*) T ), 
where T' = ^ ^fl^fJ 2 ,^^ is a constant. The FP equation that would result from the above 

<ll(Z)Xl2(/)> ^ 

Langevin dynamics would be of the form 

= ~{C.nx)P(x, t)) + T'K B T^- 2 P(x, t). (12) 

The normalized steady state solution for this FP is easily identifiable as 

P(x) = A^exp( t , KbT )■ (13) 

Thus, we identify the constant in Eq.9 which makes the exponent non-dimensional as 
C = 1/T'KbT. Note that, for a constant damping we immediately recover the standard 
Boltzmann distribution function from our general expression for P(x). 



DISCUSSION 



The first consequence of having this probability distribution is, in general, X =< 
(xF(Z) >^ (compare the forms of and (z in terms of I\s). Thus, such a system 
can have a uniform centra of mass motion in the homogeneous bath. The homogeneity 
of bath-space does allow this unlike the conformation space which is inhomogeneous due 
to the presence of the internal field. Generally, people deny the presence of a CM motion 
under equilibrium conditions by mixing up ideas of the current in the presence of a global 
symmetry breaking field which appears for a non-equilibrium steady state with associated 
entropy production. But, in homogeneous space a CM velocity will not be associated with 
any entropy production because its a system in neutral equilibrium where some initial ve- 
locity gained while equilibrating can, in general, be maintained if symmetry breaking of 
the internal forces by damping (external force) can balance the overall damping on the CM 
motion. A particular case has already been discussed in 

Let us take a standard form of Langevin dynamics for a unit mass particle (for simplicity) 
with the conformation dependent damping and generalized FD relation as 



dU(x) „, , / 2K B T . . 

The corresponding velocity distribution being P(v(x)) = N exp 2 K B TT(x)/<r(x)> which is 
similar in form to the expression we have taken in terms of (s and T f (which are basically 
inverse damping), one can do a local velocity average to get < v(x) 2 > with the above 
mentioned distribution to convert the strength of the noise to get the relation 

r{x)2 TT^ = T{x) < v{x)2 >; (15) 

which is nothing but a statement of maintenance of local balance of energy. This readily gives 
a justification of the generalization of the Maxwell distribution of the system in accordance 
with the the generalized FD relation. This also lends support to the generalization of the 
FD relation on other hand showing that it keeps the local energy balance intact. 

Let us try to understand the most striking advantage of such a situation of having gen- 
eralized distribution for conformation dependent damping. Its just evident from the very 
sight of the probability distribution that an energy barrier overcoming can be aided by a 
large local damping which advantage is absent in the standard Boltzmann distribution. May 



be such a mechanism is in use in the problem of protein folding! May be its in use in a 
whole lot of other biological and glassy systems to help overcome the energy barriers and 
the landscape is not actually as rugged as it appears from the simple view of the potential 
itself! In experiments, actually measured deviation of the local probability distributions of 
the various conformations from the standard Boltzmann one can verify this present approach 
of generalization of the theory which probably so far not been done because of over reliance 
on the FD relation with a constant damping. 

The generalization of the FD relation in the present context, which was initiated by 
getting a zero average velocity in the inhomogeneous conformation space of the system 
in order to attain an equilibrium distribution, indicates a generalization of the standard 
Maxwell-Boltzmann distributions. This is not surprising, because, any system in equilibrium 
or in other words having equilibrium fluctuations would encounter an equilibrium damping. 
So, no system in equilibrium with a heat bath is actually a Hamiltonian system in the 
microscopic sense. It appears Hamiltonian in the presence of FD relation and at the scales 
on which the FD relation is defined. Moreover, the standard derivation of the FP with a 
uniform damping constant or just straightforward replacement of T by T(x) (which does 
not ensure zero current in the present model considered) gives us the clue as to why it 
does not show up in the equilibrium probability distribution, because, it cancels out in the 
numerator and the denominator. This process gives us the Boltzmann distribution which is a 
consequence of considering a Hamiltonian flow in the phase space (consequence of Liouville's 
equation). But, a nonuniform damping, which most generally can be the case, actually calls 
for a generalization and in the present work we have captured such a situation by suitably 
generalizing the FD relation. 

Now, this generalization also imposes a problem. The damping now is not a constant 
only dependent upon the unknown degrees of freedoms that steal energy from the degrees 
of freedom are actually considered, rather, it being as well a function of the explicitly 
considered degrees of freedom there should be a prescription for finding it, otherwise, we 
lose predictability. This is where one should now look at to identify some general guiding 
principles for the conformation dependence of the damping where the unknown degrees of 
freedom part can feature as an adjustable constant. 
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